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ABSTRACT 

The search-volume corrected period distribution of contact binaries of the W UMa 
type appears to reflect primarily the constant number ratio of ~ 1/500 to the number 
of stars along the Main Sequence; there exist no evidence for angular momentum 
evolution. The maximum in contact binary numbers is located at shorter periods 
than estimated before, P ~ 0.27 d. The drop in numbers towards the cut-off at 
P ~ 0.215 — 0.22 d still suffers from the small number statistics while the cut-off 
itself remains unexplained. Only one out of seven short-period ASAS variables with 
P < 0.22 d have been retained in the sample considered here within 8 < < f 3; this 
short-period field-sky record holder at P = 0.2178 d should be studied. 
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1 INTRODUCTION 

For Main Sequence binary stars, the orbital period (P) 
distribution tells us primarily about the binary forma- 
tion processes and the angular momentum {H) distribu- 
tion among the pre-stellar fragments. These processes ap- 
pear to produce a featureless, (almost) scale- free distribution 
n(P)d P PC (l/P)dP with a ver y mild maximum at ~ 180 
years (|Duauennov. Mavoj Il99ll . =DM91). The dominance 
of the und erlying l o garith mic distribution was particularly 
stressed by iHeacoxl 1 2000l ): in fact, sometimes in theoretical 
simulations, a flat distribution n(logP)dlogP = const is 
assumed. 

Contact binaries of the W UMa type do not obey the 
flat logarithmic period distribution: Not only a very sharp 
edg«El to the distribution appears at about 0.215 - 0.22 d, 
but a strong maximum in the period distribution is present 
very close to the cut-off at very short periods. This is illus- 
trated in Figure [T] showing data for the best observed 352 
sh ort-period (P < 1 day) binary systems from the catalogue 
of iPribuUa et all |2003). The sharp edge is particularly well 
deflned in the n{P)dP distribution (left panel), but - in view 
of what has been said above - a better way to analyze the 
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1 Among field star s, CC Com with P = .2207 d has the shortest 
period. V34 in the IWeldrake et al. I ||2004| ') survey of the globular 
cluster 47 Tuc has an even shorter period of P = 0.2155 d; contact 
binaries in glob ular clusters are expected to be smaller and tighter 
llRucinskill2000h . 




Figure 1. The period distrib ution based on the c atalogue data 
for bright contact binaries in iPribuUa et al. I (120031) . It is shown 
as n{P)dP (left panel) and in n(log P) dlog P (right panel) dis- 
tributions. 



numbers may be to use logP units (right panel). Although 
the short period cut-off is then less sharp, a well defined 
maximum is still there. However, the catalogue data may 
be affected by several selection effects. First of all, the very 
steep period - luminosity relation (see below) results in a 
strong dependence of the search volume on the period. In 
addition, random discovery conditions and preferential at- 
tention of variable-star observers produce additional, totally 
uncharacterized selection effects. The period range of 6 to 
12 hours, exactly the region of interest here, is particularly 
prone to various subjective preferences. 
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An attempt was made by iRucinskil l|l992l ) to explain 
the period cut-off as a "full convection limit" for low mass 
stars. Although this gave some insight into the physics of 
the least massive contact binaries, it failed to explain the 
existence of the cut-off. Because the W UMa-type contact 
binaries are very magnetically active, an explanation for the 
cut-off may be related to their activi ty, bu t it is not clear 
how. Following the pioneering work of lVilhul (|l982h and sev- 
eral observational evidences of "saturated" magnetic (coro- 
nal, chr o mosp heric) activity at very short rotation periods, 
IStepieril l|2006l ) attempted to explain the period cut-off by 
the magnetic- wind driven angular momentum loss (AML), 
operating in such a saturated regime. He found that at the 
shortest periods, the AML rate depends primarily on the bi- 
nary moment of inertia (which can vary in contact binaries 
due to possibility of the mass transfer between components). 
This results in a progressive decay of the AML rate with the 
shortening of the period so that the period evolution takes 
progressively longer time. The period cut-off would be then 
due to the finite age of the binary population of several Gyr. 
Note that the tendency of AH/H ~ const would produce a 
logarithmically flat distribution of periods, AP/P ~ const, 
while a decreasing AH/H would lead to a pile-up at short 
periods. 

In this paper, we analyze the shape of the pe- 
riod distribution close to the period cut-off. To avoid 
the most obvious selection effects, the uniform and well 
characterized ASAS s ample l|Poimahski. MacieiewskJ 



2005 



2ooe 



iPaczvhski et al.ll2006l) is used, the same as in Rucinski 
Paper I) where it was used to estimate the local fractional 
number density of W UMa systems at ~ 1/500 of the FGK 
dwarfs. The same approach was utilized there so that the 
reader is referred to previous paper for details. 



2 THE ASAS SAMPLE 

The material for this study is the ASAS sample discussed 
in Paper I, consisting of 5381 binaries and apparently uni- 
form in terms of detection selection effects within the range 
8 < y < 13. Although many small amplitude binaries are 
included in the sample, down to amplitudes of ~ 0.05 mag, 
the amplitude detection threshold for completeness is rela- 
tively high for this sample at about 0.4 mag, necessitating 
a correction for missed low-amplitude binaries of approxi- 
mately (3 ± 1) times (Paper I). 

The ASAS photometric data are in only one V band. 
For that reason the My = M y (log P, B — V) calibration of 
iRucinski fc Duerbeckl (|l997h could not be used for distance 
estimates. This forced utilization of a sub-sample of 3374 
shortest-period binaries where a simplified period - absolute 
magnitude calibration appears to be valid (Paper I). This 
way the dominant selection effect of the search volume being 
directly dependent on the binary period could be taken into 
account. The period dependence in the calibration is very 
steep, Mv ~ —1.5 — 12 x logP; it certainly does not apply 
for log P > —0.25 (or P > 0.562 d) where the luminosity is 
seen to depend on both, the period and the colour index. 

The selection of the ASAS sample used in Paper I was 
not altered in any way except for a particular attention to 
existence of contact binaries at P < 0.22 d. The on-line 
version of the ASAS database was carefully scrutinized (Ap- 
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Figure 2. The volume-corrected period distribution based on the 
ASAS data, taking into account the distance selection effects. The 
error bars are derived from the Poisson statistics uncertainties 
based on the number of objects, as given in Table [T] The broken 
line gives the DM91 distribution scaled with the assumption of 
the relative frequency of occurrence of one W UMa binary per 
500 FGK dwarfs. 



Table 1. The volume-corrected period distribution Dp derived 
from the ASAS sample and expressed as number of contact sys- 
tems per interval AlogP(d) = 0.025 in units of 10~^ stars pc"''. 
n is the number of binaries per interval; it can be used to derive 
the Poisson uncertainty of Dp. The first column gives the low 
edge of each bin. 



log 


p 


ri (pc) 


r2 (pc) 


n 


Dp 


-0 


675 


21.88 


190.55 


3 


0.423 


-0 


650 


25.12 


218.77 


8 


0.746 


-0 


625 


28.84 


251.19 


12 


0.739 


-0 


600 


33.11 


288.40 


38 


1.546 


-0 


575 


38.02 


331.13 


74 


1.989 


-0 


550 


43.65 


380.19 


92 


1.634 


-0 


525 


50.12 


436.52 


130 


1.525 


-0 


500 


57.54 


501.19 


155 


1.202 


-0 


475 


66.07 


575.44 


233 


1.193 


-0 


450 


75.86 


660.69 


260 


0.880 


-0 


425 


87.10 


758.58 


279 


0.624 


-0 


400 


100.00 


870.96 


312 


0.461 


-0 


375 


114.82 


1000.00 


247 


0.241 


-0 


350 


131.83 


1148.15 


234 


0.151 


-0 


325 


151.36 


1318.26 


184 


0.078 


-0 


300 


173.78 


1513.56 


161 


0.045 


-0 


275 


199.53 


1737.80 


167 


0.031 



pendix|X]| and - indeed - from 7 systems previously included 
in Paper I, only one has remained as a genuine contact bi- 
nary (083128-1-1953.1, P = 0.217811 d); this is currently the 
shortest period record holder among galactic field W UMa 
binaries. 



3 THE MAIN RESULTS 

In Paper I the period distribution for all contact binaries 
of the ASAS survey was shown (Fig. 3 in that paper). It 
contained the strong period-distance selection effect which 
we now attempt to eliminate. For each binary of known P 
and V, the simplified calibration My = —1.5 — 12 x logP 
was used to derive My and then the distance r (in pc) from 
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Figure 3. The volume-corrected period distribution based on 
the ASAS data, taking into account the distance selection effects. 
This is the same diagram as in Figure (2] but in linear units of 
the period. The error bars are derived from the Poisson statistics 
uncertainties, as in Table |2] Note that there are only 3 systems 
in the first bin of 0.200 < P < 0.225 d. 



Table 2. The volume-corrected period distribution D'p derived 
from the ASAS sample and expressed as number of contact sys- 
tems per interval AP(d) = 0.025 in units of 10~® stars pc~^. n 
is the number of binaries per interval. The first column gives the 
low edge of each bin. 
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ri (pc) 


'■2 (pc) 


n 




0.200 


22.14 


166.91 


3 


0.630 


0.225 


28.51 


221.43 


18 


1.619 


0.250 


35.84 


285.14 


77 


3.243 


0.275 


44.17 


358.43 


129 


2.735 


0.300 


53.52 


441.66 


179 


2.028 


0.325 


63.94 


535.20 


257 


1.636 


0.350 


75.45 


639.39 


304 


1.135 


0.375 


88.09 


754.53 


314 


0.713 


0.400 


101.89 


880.94 


321 


0.458 


0.425 


116.87 


1018.90 


245 


0.226 


0.450 


133.07 


1168.72 


207 


0.127 


0.475 


150.50 


1330.65 


173 


0.072 


0.500 


169.19 


1504.97 


131 


0.038 


0.525 


189.18 


1691.93 


129 


0.026 


0.550 


210.48 


1891.78 


74 


0.011 



V — Mv = 5 logr — 5. Then the numbers of binaries, n, 
in intervals of A log P were counted within the limits of 
8 < y < 13 corresponding to the respective distance lim- 
its ri and r2 (Table [T|. By dividing the numbers n by the 
volume of the corresponding shells defined by ri and r2, 
and taking into account the sky coverage of the ASAS sam- 
ple (73.5%) and the amplitude-related under-counting by a 
factor of 3x, we could derive the volume-corrected period 
distribution Dp. (This quantity was called a "period func - 
tion" in the analysis of the OGLE sample in lRucinskil (|l998l ). 
see Fig. 3 there). Dp is expressed in numbers of stars per 
counting interval AlogP and per unit of volume in pc^. 
This is basically the same approach at that used in Paper I 
for derivation of the luminosity function. Even Figure 4 of 
Paper I can be of interest here if its horizontal axis with 
My units is converted into the logP units (note for exam- 
ple the vertical strip of missing systems with P ~ 0.5 d at 
My — +2.1, an effect most likely caused by the sparse and 
somewhat regular data taking of typically 2-5 observa- 
tions per night and by associated difficulties with detection 
of periods commensurable with one day). The results for 
the counting interval of AlogP(d) = 0.025 are shown in 
Figure [1 

The period distribution in Figure [2] shows a very strong 
peak at about logP ~ -0.57 or P ~ 0.27 d. Thus, the 
maximum appears at periods shorter than for the catalogue 
data (Figure [l]) which directly shows the systematic effects 
of the sampling volume being dependent on the period. One 
also sees a drop in numbers shortward of the maximum, 
but the Poisson uncertainty there is large. Even with the 
search limits extending beyond V — 13, the ASAS survey 
has small number of systems with P 0.27 d. An artificial 
gap at 0.25 d, similar to that definitely present at 0.5 d, may 
additionally depress the distribution at very short periods. 

4 COMPARISON WITH THE DETACHED 
MAIN SEQUENCE BINARIES 

The volume-corrected period distribution for contact bina- 
ries derived above. Dp, can be compared with an extension 



of the Main Sequence (MS) period distribution for detached 
binary stars of DM91. This involves a rather extensive and 
perhaps risky extrapolation into the domain of P < Id 
where the DM91 survey contains only one object. 

The MS binary distribution has the shape Dms oc 
exp(-(logP - logPo)^/(2cr|>)) with logPo = 4.8 and crp = 
2.3, where the period unit is one day. After normalization of 
this distribution to the integral of unity, the distribution was 
multiplied by 0.53 to represent the binary fraction among 
MS stars (DM91). Such a normalized distribution must be 
further adjusted for a direct comparison with the contact- 
binary Dp distribution. Integration of Dp (as in Fig. [2} over 
the whole range of periods, after accounting for the contact- 
binary occurrence fraction of 1/500 gives the total number 
density of the MS stars in the ASAS survey volume con- 
sidered here of 6750 stars/ (lO^pc^); this was the multiplier 
used to derive the line shown in Figure [2] to represent the 
expected MS binary frequency. An entirely independent esti- 
mate of the multipli er based on the local luminosity function 
(jWielen et al.lll983l l gives 9350 stars/ ( 10 V^)- The ~ 1/3 
discrepancy in these numbers is a fair representation of the 
uncertainty inherent here. 

As discussed in Section [1] the expected Main Sequence 
period distribution Dms, as shown in Figure [2] is almost 
flat within the contact binary range, ft is also drastically 
different from what we see for contact binaries. Apparently, 
contact binaries are very common and dominate in numbers 
in the period range of 0.25 < P < 0.35 d, but are relatively 
infrequent for longer periods. Presumably, other binaries are 
common there, but any quantitative studies of this would 
have to address the very different discovery selection effects 
for contact, semi-detached and detached short-period bina- 
ries. We note that ILucvI (j 19761 ) bravely attempted to use 
the General Catalogue of Variable Stars to address exactly 
the same point, but his data were very poor at that time. 
Surveys such as ASAS or many planned ones can be used 
to study the metamorphoses of short-period binaries as they 
evolve into contact; here we only point interesting potential 
for research in this area. 

In Paper I it was found that the luminosity function for 
contact binaries scales as 1/500 of the MS function for the 
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range 1.5 < My < 5.5 (roughly A3V to G6V). Thus, the 
numbers of contact binaries simply vary along the MS as 
the numbers of stars increase towards low luminosities; the 
short-period cut-off is located within the 5.5 < Mv < 6.5 
bin of G6V to K2V stars where the number fraction drops 
from 1/500 to approximately 1/1000. The continuously vary- 
ing mass with the period for contact binaries is a very dif- 
ferent situation than for detached binaries where DM91 as- 
sumed the same masses and thus where the period distri- 
bution could reflect the angular momentum content and its 
changes. The flat shape of the n(log P)ci log P distribution 
is therefore not necessarily expected to apply to contact bi- 
naries and indeed is not observed. (For completeness, Fig- 
ure[3]and Table[2]give the period distribution in linear units, 
D'p ). Since the numbers of contact binaries reflect mostly 
the mass distribution, not much can be said about the evo- 
lution of the angular momentum. It is al so hard to pr ove 
or disprove the saturated AML evolution (|St(?pierill2006l ) on 
the basis of the period distribution. 



5 CONCLUSIONS 

We found that the search-volume corrected period distribu- 
tion for W UMa binaries is a direct result of their luminosity 
function following the Main Sequence luminosity function 
over the range of 1.5 < My < 5.5. The maximum of the pe- 
riod distribution occurs at periods close to 0.27 d {Mv — 6), 
at the point where the W UMa luminosity function has a 
point of inflection and stops following the MS numbers. The 
maximum is located at shorter periods than estimated pre - 
viously from the catalogue data (|LucvII197^ : iRucins kil ll992l ) 
at about 0.35 d; see also Section [U The new determination 
based on the ASAS survey and based on 3374 objects is 
clearly superior to the one based on 98 objects of the OGLE 
survey (|r ucms kil ll998l l. The OGLE sample analysis appears 
to have suffered from an unaccounted image blending which 
led to erroneous depth limits and, in consequence, to an 
exaggerated frequency of occurrence of W UMa binaries of 
~ 1/130 in place of ~ 1/500. But the shape of the period 
distribution was probably unaffected. 

The statistics at shorter periods beyond the maximum 
is still very poor. Apparently, even a survey such as ASAS 
with the coverage of almost 3/4 of the sky and extending 
to 1/ ~ 13 still has too few objects to establish the statis- 
tics at short periods; deeper all sky surveys are needed. 
Nevertheless the ASAS data have improved our knowl- 
edge of the short-period end of the distribution and of its 
maximum. The shortest-period known fleld system, ASAS 
083128+1953.1, with P = 0.2178 d needs investigation; 
other ASAS variables with P < 0.22 d which were previously 
considered in this context are either not contact binaries or 
have poor light curves at their brightness levels beyond the 
V = 13 limit. 

The use of the n(logP)dlogP or n{P)dP forms of the 
period distribution does not seem to make such an impor- 
tant distinction for contact binaries as for Main Sequence 
stars of same mass. The strong maximum in numbers ap- 
pears at short periods in both of these forms and is mostly 
related to the stellar mass distribution along the Main Se- 
quence. It should be noted that - unlike detached binaries 
- contact binaries can respond to the angular momentum 



Table 3. Systems with periods < 0.22 d. 



ASAS 


P(d) 


V 


B - V 


Comment 


042606+0126.5 


0.148860 


12.4 


0.27 


S Set, 0.074430 d 


071829-0336.7 


0.211249 


13.8 


? 


> 13 mag, poor I.e. 


083128+1953.1 


0.217811 


10.3 


1.03 


genuine W UMa 


113031-0101.9 


0.213135 


13.4 


? 


> 13 mag, poor I.e. 


162155-2128.7 


0.209918 


10.1 


0.47 


too blue, ampl. 0.06 


174930-3355.4 


0.186629 


13.8 


? 


> 13 mag, poor I.e., new V 


201354-4633.4 


0.147786 


9.6 


0.17 


too blue, ampl. 0.08 



loss by an internal mass transfer towards more dissimilar 
masses (i.e. a decrease in the mass ratio) with the associ- 
ated lengthening of the period. However, this is not directly 
visible in the ASAS period distribution at the present level 
of the statistical accuracy. 

The author would like to express his thanks to Dr. G. 
Pojmahski for making ASAS such an useful variable-star 
study tool and to Dr. K. St§pieh for reading a draft of the 
paper and for very useful suggestions. Thanks are also due 
to the reviewer for succinct and focussed suggestions. Sup- 
port from the Natural Sciences and Engineering Council of 
Canada is acknowledged with gratitude. 
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APPENDIX A: CONTACT BINARIES WITH 
PERIODS SHORTER THAN 0.22 DAY IN THE 
ASAS SURVEY 

The seven systems of the ASAS sample with P < 0.22 d 
(Table|3]l were checked if they are genuinely contact binaries. 
The on-line version of the ASAS-3 catalogue was consulted 
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for updates and the Tycho-2 catalogue l|H0g et al.ll200Cf ) was 
used as a source of B ~ V colour index data. 

It appears that the system 042606+0126.5 has been re- 
classified in the ASAS catalogue to a 5 Set pulsating star 
with the period 2 times shorter than estimated before. Also, 
the ASAS tables give V = 12.6 for 174930-3355.4 while its 
rather poor light curve clearly shows V = 13.8; thus, the star 
is fainter than the adopted completeness limit of V — 13. 
Two more stars are fainter than this limit and also have poor 
light curves, 071829-0336.7 and 113031-0101.9. 

Short period contact binaries must be red because they 
consist of mid-K dwarfs. Thus, a colour index which is not 
sufficiently red {B~V > 0.8) can be used to reject S Set and 
j3 Cep pulsating stars. This criterion was used to eliminate 
162155-2128.7 and 201354-4633.4. 



